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We show that the presence of a p ~ r~ 3 ^ 2 dark matter overdensity can be robustly predicted at 
the center of any galaxy old enough to have grown a power-law density cusp in the stars via the 
Bahcall-Wolf mechanism. Using both Fokker-Planck and direct iV-body integrations, we demon- 
strate collisional generation of these dark matter "crests" (Collisionally REgenerated STtructures) 
even in the extreme case that the density of both stars and dark matter were previously lowered 
by slingshot ejection from a binary supermassive black hole. The time scale for collisional growth 
of the crest is approximately the two-body relaxation time as defined by the stars, which is < f 
Gyr at the centers of stellar spheroids with luminosities L < 10 9 5 Lq, including the bulge of the 
Milky Way. The presence of crests can robustly be predicted in such galaxies, unlike the steeper 
enhancements, called "spikes", produced by the adiabatic growth of black holes. We discuss special 
cases where the prospects for detecting dark matter annihilations from the centers of galaxy haloes 
are significantly affected by the formation of crests. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

While the evidence for a dynamically significant com- 
ponent of dark matter (DM) on cosmological scales is 
compelling, the nature of the DM is unknown, and its 
distribution on sub-galactic scales remains uncertain. A 
widely discussed DM candidate is the supersymmetric 
neutralino 0,0; the presence of which might be detected 
indirectly through the products (gamma rays, neutrinos, 
anti- matter) of its self-annihilations p|. |_l|| . In the case of 
photons, the annihilation signal is simply proportional to 
the square of the DM density p x integrated along the line 
of sight, and most discussions of indirect detection have 
focussed on the centers of galaxies, including the Milky 
Way galaxy, where the DM density is likely to be highest 

Supermassive black holes (SBHs) are believed to be 
generic components of galactic nuclei |Toj and are ex- 
pected to strongly influence the distribution of mass 
(stars, DM) at distances < 77, from the SBH 11], where 
rh is the gravitational influence radius, defined as the ra- 
dius within which the gravitational force from the SBH 
dominates that from the stars. In the case of the Milky 
Way SBH, m k, 3 pc. In one widely discussed model |12| . 
a so-called spike forms around the SBH at r < as it 
grows adiabatically. The density in such a spike is a steep 
function of radius near the SBH, p^ ~ r -7 , 7 > 2 imply- 
ing a large DM annihilation rate |l2Lll3l[l^ . However the 
formation of such spikes requires finely-tuned initial con- 
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ditions and the spikes are easily destroyed 0, 0] . 
Furthermore, the stars would react in the same was as the 
DM to spherically-symmetric growth of a SBH ^lj , and 
in the galaxies with long central relaxation times where 
steep stellar cusps could persist for 10 Gyr or longer, none 
are seen ^lj. In such galaxies, stellar density profiles are 
typically flat at r < 77,, believed to be a consequence 
of the "scouring" effect of binary SBHs during galaxy 
mergers 0] . The density of DM at the centers of these 
galaxies would also presumably be low. 

Steeply-rising stellar densities are instead observed at 
r % r h in the bulge of the Milky Way and possibly in 
M32, a nearby dwarf elliptical galaxy pj||2(J . In both of 
these dense, compact stellar systems, the two-body relax- 
ation time near the SBH - the time for stars to exchange 
orbital energy via gravitational encounters - is < 10 10 
yr, short enough for the formation of a Bahcall-Wolf plf 
("collisional") density cusp in the stars around the SBH, 

~ r~ 7 / 4 . In dense galaxies like these, a Bahcall-Wolf 
cusp in the stars can even re-form after being destroyed 
by a binary SBH [12. 

In this paper, we discuss the evolution of the DM den- 
sity in a nucleus that grows a collisional cusp in the 
stars via the Bahcall-Wolf mechanism. The DM parti- 
cles are essentially collisionless, but they scatter off of 
stars [23l |2J, |25j , forming a p x ~ r -3 / 2 density "crest" 
(Collisionally REgenerated STructure) near the SBH in 
roughly one stellar relaxation time. Remarkably, as we 
show, this is true even in the case that the galaxy core 
was previously "scoured" by a binary SBH: DM particles 
are scattered by stars into regions of phase space corre- 
sponding to tightly-bound orbits around the SBH that 
were previously depleted by the binary. 

In Sects. II and III we present Fokker-Planck, as well 
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as direct TV-body, integrations of a combined, star+DM 
system around a central point mass that demonstrate the 
formation of crests on roughly a star-star relaxation time 
scale. Sect. IV discusses the envirnomental conditions 
necessary for the formation of crests; we show that these 
conditions are likely to be satisfied in stellar spheroids 
comparable in luminosity to that of the Milky Way or 
fainter, allowing us to robustly predict the presence of 
DM crests in these systems. In Sect. V we discuss spe- 
cial cases where the prospects for detecting dark matter 
annihilations from the centers of galaxy haloes are signif- 
icantly affected by the formation of crests. Our conclu- 
sions are summarized in Sect. VI. 



II. FOKKER-PLANCK TREATMENT 

An approximate description of the combined evolution 
of stars and DM at the center of a galaxy is provided 
by the isotropic, multi-mass Fokker-Planck equation [2^, 
123,123. Let f(E, to, t)dm be the number density in phase 
space of objects (stars, DM particles) in the mass range 
to to to + dm; E = —v 2 /2 + <f> > is the binding energy 
per unit mass and $ = — <p is the gravitational potential, 
assumed fixed in time. (In what follows, changes in / are 
only significant within the SBH's sphere of influence and 
the assumption of a fixed potential is reasonable.) Then 
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e.g. g is the phase-space mass density. 

Consider now a nucleus containing just two compo- 
nents, stars of mass and DM particles with an un- 
specified range of masses such that m x <C to*. Assume 
further that g x <C g* ■ Taking the first moment of Eq. i|la|) 



over mass then yields the two evolution equations 
1 d 
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with diffusion coefficients 
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The DM particles, being of negligibly small mass, do not 
self-interact gravitationally and they evolve solely due 
to heating by (i.e. scattering off of) the stars |23L l24j. 
As a result, the characteristic time for change in cither 
<?* or g x is the same, (47iTm*g*) , equal to within a 
constant factor to the standard two-body relaxation time 
T r defined by the stars alone p7|: 
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with v rrns the mean square velocity. In a time ^ E ri the 
stellar distribution approaches its steady-state form near 
the SBH, 



g^E) ~ E x l\ p*(r)~r- 
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the so-called Bahcall-Wolf [21| solution. The steady-state 
solution for the DM is obtained by setting dg x /dE = 0, 
or 



9x( E ) ~ const., p x (r) - r 



-3/2 
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|24l I25I l29j . Both of these solutions assume a non- 
evolving phase space density far from the SBH; in re- 
ality, the DM density will drop after the formation of 
the crest, due to ongoing heating by the stars |2^, and 
under certain circumstances the stellar density may con- 
tinue to evolve as well (* ne functional form Q of g x 
is unaffected by such evolution and we ignore it in what 
follows). While the solution J7J for g x at r <C is in- 
dependent of , the time required to attain this density 
profile, and the subsequent rate of heating by the stars, 
do depend on g*. 

Fig. n shows time-dependent solutions to Eqs. 13al) - 
(|4b|) . Initial density profiles for both the stars and DM 
were assigned as in Ref. [30|: these models have p oc r 1 
near the SBH, and we chose 7 = 0.5, the flattest density 
profile consistent with an isotropic phase-space distribu- 
tion in a 1/r potential. (The density at large radius in 
this model falls off more steeply than the p x ~ r^ 1 depen- 
dence predicted in standard DM halo models; however we 
are concerned here with the evolution only at very small 
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FIG. 1: Solutions of the Fokker-Planck equations 13aH - 14b|l that describe the joint evolution of stars and dark matter around a 
black hole due to star-star and star-DM gravitational encounters. Length unit rn is the radius containing a mass in stars equal 
to twice the black hole mass at t — (roughly 3 pc at the Galactic center). Density is in units of its initial value at r^. Curves 
show the stellar (left) and dark matter (right) density profiles at times (0, 0.2, 0.4, 0.6, 0.8, 1.0) in units of the initial relaxation 
time (Eq. at r = ry t . Dashed lines are the "steady-state" solutions, Eqs. © and (J. 
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FIG. 2: Evolution of the stellar (solid lines) and dark matter 
(dashed line) densities at radii of (10~~ 5 , 10~ 4 , 10 _3 )r; l in the 
Fokker-Planck integration of Fig. 1. Densities are normalized 
to their values at t — 0. 



radii.) The unit of time in Fig. Q] is the relaxation time 
for the stars at the influence radius r^, defined in the 
standard way as the radius containing a mass in stars 
equal to twice M». In principle, r/, so defined changes 
with time in response to the changing stellar density, but 
this effect is almost negligible and we ignore it in what 
follows. 

In a time ~ 0.5T r (rh), the stars are seen to attain 



the Bahcall-Wolf profile at r < 0.2r^. The evolution of 
the DM is more complex. A p x ~ r~ 3 / 2 crest inside 
r ~ O.lrh is formed in roughly the same time. However 
the amplitude of the crest drops thereafter as the DM is 
heated by the stars (Figs. : by a factor 1/e in a time 
of ~ \.2T r (rh) after peak density, and by a factor 1/e 2 
in a time of ~ A.hT r (rh)- Evolution of the DM density 
profile after crest formation is approximately self-similar, 
p x (r,t) « p xfi {r)G{t/T r ) at r < r h , with dG/dt < 0. 



III. iV-BODY TREATMENT 

The Fokker-Planck equations (|3afl - i|4b(l embody a num- 
ber of approximations (isotropy, small-angle scattering, 
uncorrelated encounters, fixed gravitational potential, 
etc.) that may be violated in real stellar systems. Fur- 
thermore the initial conditions of Fig. ^ are ad hoc. We 
therefore carried out a direct TV-body integration of stars 
and massless particles around a point mass, starting from 
initial conditions that realistically represent the center of 
a gal axy after a binary SBH has created a low-density 
core 0]. 

The iV-body code was adapted from </?GRAPE [3l| . 
a direct-summation code that advances particles via a 
fourth-order integrator and computes gravitational forces 
via calls to special-purpose accelerator boards called 
GRAPEs [l^]. The code was modified to run in serial 
mode using a single GRAPE-6 computer which has an 
onboard memory limit of ~ 256K particles and a speed of 
~ 1 Tflops. The code was also modified to include mass- 
less (DM) particles; since these do not "see" each other 
gravitationally, they need not all be loaded into GRAPE 
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FIG. 3: Direct JV-body integration of "star" (left) and "dark matter" (right) particles around a single massi ve p article ("black 
hole"), starting from a galaxy model in which a pre-existing binary black hole had created a low-density core |22j. Thin, normal 
and thick curves show the density profiles at times (0, 0.25, 1.0) in units of T r (r^). Dashed lines are as in Fig. 1. 



memory simultaneously, allowing the GRAPE's memory 
limit to be circumvented. We used AT* = 1.2 x 10 5 "star" 
particles and N x = 2.4 x 10 6 "DM" particles. 

Initial conditions for the TV-body integration were 
adapted from Run 8 of Ref. In that paper, JV- 

body simulations were used to follow the formation of 
low-density cores via inspiral of a SBH into the center 
of a galaxy containing a second, more massive SBH; in 
Run 8, the binary mass ratio was 1 : 4. We replaced the 
two massive particles at the final time step of Run 8 by 
a single particle (the SBH) having their combined mass 
(0.0125 in units of the total galaxy mass), and positioned 
this particle at the center of mass of the pre-existing bi- 
nary. Massless (DM) particles were then added via a 
bootstrap algorithm: a star particle was chosen at ran- 
dom; a DM particle was placed randomly on a sphere 
with radius equal to that of the star particle; the DM 
particle was assigned a velocity with random direction 
subject to the constraint that the radial and tangential 
components were equal to those of the star particle. This 
scheme produced a model in which the DM had the same 
initial phase-space distribution as the stars; in particu- 
lar, both components had essentially flat central density 
profiles (Fig. |3J). In addition, the velocity distributions 
at r < rh were biased toward circular motions, a result 
of gravitational slingshot ejection of particles on radial 
orbits by the binary. A-body integration of this model 
until a time ~ T r (rh) using (^GRAPE required ~ 43 d 
on the GRAPE-6 special-purpose computer, for a total 
of ~ 6.6 x 10 17 floating-point operations. 

Fig. |3| shows the results. The star particles form a 
Bahcall-Wolf cusp in a time ~ 0.5T r (r/j). The response of 
the DM particles is also consistent with what was found 
in the Fokker-Planck integration (Fig. ^| : formation of a 
crest at r < O.lr^ and a gradual drop in p x as the DM 



is heated by the stars. The DM density slope at r < 
O.Olr/j is not well constrained due to the finite number of 
particles but is consistent with p x oc r~ 3 / 2 . Henceforth 
we will assume that the Fokker-Planck solution correctly 
describes the radial dependence of the DM density at 
radii too small to be resolved via the A-body integration. 

Both the size of the core formed by an inspiralling 
SBH (measured in units of the binary's mass) and the 
time scale to regrow a stellar cusp (measured in units of 
Tri rpX ] are almost independent of the binary mass ratio 
|18l l22l ] and so the results of our single A-body integra- 
tion should be representative of the majority of galactic 
nuclei that formed via mergers, if the length and mass 
are scaled to and M, respectively. 



IV. GALAXY PROPERTIES RELEVANT TO 
THE FORMATION OF CRESTS 

Several basic conditions must be satisfied for a DM 
crest to form in a galactic nucleus. (1) The two-body 
(star-star) relaxation time must be short. (2) The nu- 
cleus must contain a massive BH. (3) For the annihila- 
tion signal to be detectable, there must be a significant 
amount of DM on scales < r/,. In this section, we re- 
view current knowledge concerning physical conditions 
in galactic nuclei and discuss which kinds of galaxies are 
most likely to harbor DM crests. 

A. Relaxation times 

The two-component models presented above demon- 
strate that a stellar cusp and its associated DM crest 
are generated in a time of ~ 0.5T r (r/ l ). (The single- 
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FIG. 4: Relaxation times at the SBH's influence radius 
in Virgo cluster galaxies (circles) and the Milky Way (star) 
vs. the total blue luminosity of the galaxy (in the case of the 
Milky Way, the bulge). A stellar mass of IMq was assumed 
when computing T r . Filled circles: Virgo cluster galaxies in 
which rh is resolved, i.e. rn > 0.1". Open circles: Virgo 
cluster galaxies in which is unresolved. Stellar luminos- 
ity profiles were taken from Ref. [33j| and black hole masses 
were computed from the M. — a relation [To||. except in the 
case of the Milky Way for which M, has been directly deter- 
mined from stellar orbits |34|. Horizontal dotted lines are at 
T r (r h ) = (1,2,4) x 10 10 yr. Dashed line is a regression fit to 
the data, Eq. [H] 



component iV-body studies in Refs. |22l l35l l3q reach 
similar conclusions about the time scale for stellar cusp 
formation.) Fig. 0]shows estimates of T r {rh) vs. spheroid 
luminosity in a complete sample of early-ty pe ( elliptical 
and lenticular) galaxies in the Virgo cluster |33l l37| and 
the bulge of the Milky Way. (The spheroid is defined as 
the entire luminous galaxy in the case of elliptical and 
lenticular galaxies, and as the bulge component in the 
case of spiral galaxies.) There is a well-defined trend of 
T r (vh) with spheroid luminosity L: 



T r (r h ) « 3.8 x 10 9 yr L% 



where L 9 = L/10 9 L Q . Spheroids with L < 2 x 10 9 L© 
have T r (rh) < 3 x 10 10 yr, which is short enough for 
the formation of a crest, assuming that the SBH and 
spheroid were in place at least ~ 10 10 yr ago. The Milky 
Way (marginally) satisfies this condition, and the Galac- 
tic center stellar cluster in fact has a density profile that 
is consistent the Bahcall-Wolf form [l9l l2Cj . 

Conditions for the formation of crests are relaxed some- 
what if there is a top-heavy spectrum of stellar masses 
since the DM scattering time scales as m" 1 where fa* = 
( TO *)/( m *) 24]. The stellar cusp can also evolve more 
quickly in this case [38| . 



Unfortunately, in the luminosity range most relevant 
to DM crest formation, the Milky Way bulge is the only 
spheroid near enough for a Bahcall-Wolf cusp to be de- 
tected even if present. Spheroids fainter than L sw lO 9 i0 
and outside the Local Group are always unresolved on 
scales r <Th (e.g. Fig.0}. However there are indications 
that nuclear structure begins to change systematically as 
L drops below ~ 10 9 Lq, since an increasingly large frac- 
tion of spheroids exhbit compact stellar nuclei [U • 
Only one compact nucleus is spatially well resolved, in 
the Local Group galaxy NGC 205 (L w 10 8 - 4 L Q ), and the 
relaxation time is found to drop to ~ 10 8 yr at the center 
[ioj: this value is derived assuming no SBH since NGC 
205 shows no evidence of rising velocities near the cen- 
ter [iof. Whether the compact nuclei in faint spheroids 
are consequences of their short nuclear relaxation times 
pd"! , the absence of SBHs in these galaxies 0! E3 > gas- 
dynamical processes [44|, or some combination of these 
factors is currently unclear. However the form of a DM 
crest is essentially independent of p*(r) (§2) and if the 
compact nuclei satisfy the other two conditions identi- 
fied above (presence of a SBH and a significant amount 
of DM; see discussion below) they would be expected to 
contain DM crests as well. 

If a nucleus is much older than T r (rh), a DM crest 
will form then decay in amplitude as the DM particles 
attempt to reach equipartition with the much heavier 
stars (Fig. m 24]). Fig. |2] suggests that the crest den- 
sity drops by a factor 10 from its peak value - corre- 
sponding to two orders of magnitude in the annihilation 
signal - in a time ~ 8 — 9T r (r^). The trends in Fig. 4 
are too poorly defined at low values of L/Lq to allow 
us to state clearly for which galaxies this would occur, 
but a straightforward extrapolation of Eq. JHJ) suggests 
that T r {r h ) falls below - 10 9 yr at L w 1O 87 L . If 
this inference is correct, it follows that crests would be 
present, with reasonable amplitudes, only for a fairly nar- 
row range of spheroid properties: i.e. spheroids older 
than ~one relaxation time but younger than many re- 
laxation times. Taking into account the scatter in Fig. 01 
the corresponding range in spheroid luminosities might 
be approximately 3 x 1O 8 L <L<3x 10 9 L Q . 

If low-luminosity spheroids do not contain massive 
BHs, nuclear relaxation times might generically be as 
short as those observed in NGC 205 and M33, i.e. 
10 7 — 10 8 yr. In this case, the stars would undergo core 
collapse in < 10 10 yr producing a r~ 2 - 25 density profile, 
and the DM density would be expected to evolve only 
slightly H|. 

At the other extreme of luminosity, galaxies with L > 
1O 1O L0 are always observed to have low-density cores 
with radii ~ |46l 1471 148| . and high central velocity dis- 
persions, implying long nuclear relaxation times (Fig.0J. 
In principle, these galaxies could still harbor collision- 
less DM spikes around their SBHs 0, but this seems 
unlikely: the nuclear structure of these galaxies is consis- 
tent with no dynamical evolution having occurred since 
the most recent merger "carved out" the luminous core 
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|24j . and the same merger event that created the stel- 
lar core would have destroyed a density spike in the DM 
[lfil | . Below we discuss the possibility that /ow-luminosity 
spheroids might partially retain their steep, collisionlcss 
spikes; these galaxies are unlikely to have experienced 
significant mergers and scattering off of stars need not 
completely convert such a spike into a weaker, Bahcall- 
Wolf crest in a Hubble time [24| . 

B. Black holes 

The presence of a massive BH is another necessary con- 
dition for the formation of a DM crest. Robust detection 
of SBHs is possible only in galaxies near enough that the 
stellar or gas kinematics can be resolved on scales <C r^. 
Reliable, dynamical SBH masses have been determined 
in a handful of such galaxies and the large-scale 

properties (mass, luminosity, velocity dispersion) of their 
host spheroids are found to obey tight scaling relations 
with M. [UlilEll, e.g. 

Af. i8 = (1.66±0.24)a 2 4 8 6±0 - 43 (9) 

[To|. where M. )8 = Af./1O 8 M and cr 20 o is the ID stellar 
velocity dispersion in units of 200 km s . The existence 
of these tight relations suggests that SBHs are ubiquitous 
in bright spheroids, but the faintest galaxy for which M, 
has been robustly determined is M32 (L«4x lO 8 ^©; in 
fact this galaxy is believed to be the remnant core of a 
once much brighter galaxy) and it is dangerous to assume 
that relations like © apply to fainter systems, including 
the spheroids with L < 10 9 L Q that are most relevant to 
formation of crests (Fig. |3J . 

Evidence for massive BHs in spheroids fainter than 
~ 10 9 L© - dwarf elliptical galaxies and the bulges of late- 
type spiral galaxies - comes almost entirely from the sub- 
set of nuclei that are "active," i.e. that emit a significant 
fraction of their energy non-thermally M. in these 
galaxies is determined indirectly by applying empirical 
relations established in more luminous active galaxies, 
e.g. between the size of the so-called broad emission line 
region and the nuclear continuum luminosity |5^ |. The 
presence of BHs with masses as low as ~ 10 5 M Q has 
been inferred l54|; with one exception (POX 52, a dwarf 
elliptical galaxy [53), all of the host spheroids are spiral- 
galaxy bulges. The BH masses in these galaxies appear 
to be consistent with relations like © 56j]. 

Whether BHs with M. < 10 6 M Q , are present in all 
low-luminosity spheroids is still unclear. In the Lo- 
cal Group, neither the dwarf elliptical NGC 205 (L « 
10 84 L Q ) or M33 (an apparently spheroid- less spiral 
galaxy) appear to have central BHs although the upper 
limits on M, in these galaxies (based on stellar kine- 
matics) are only marginally inconsistent with the scal- 
ing relations established in brighter galaxies [iol l57| . A 
significant fraction (~ 40%) of nearby galaxies show ev- 
idence of low-level nuclear activity but the activity need 
not be driven by BH accretion in every case |52| . Even if 



small BHs were present at one time in all low-luminosity 
spheroids, a number of processes are capable of ejecting 
them from such environments |58j . 

Some nuclei might contain binary or multiple BHs, if 
the BHs that were deposited there during a galaxy merger 
failed to coalesce. In the presence of a binary SBH, 
growth of a collisional cusp would be inhibited as the 
binary continued to eject stars and DM particles from 
the nucleus via the gravitational slingshot [5^. Time 
scales for binary SBH coalescence due to gravitational 
wave emission are longer than 10 10 yr unless the binary 
separation a drops below 

a GW ^2x 10-3 pc — i— M x 3 2i6 (10) 

where q — M<z/M\ < 1 is the binary mass ratio and 
Mi2, 6 = (Mi+M 2 )/10 6 M Q [gjj. This separation is 1-2 
orders of magnitude less than the separation at which 
the two BHs become gravitationally bound, and it is pos- 
sible for the binary to stall at a sa 3> acw (the "final 
parsec problem"). Stalling is least likely in nuclei with 
short relaxation times, since stars will scatter into the 
binary's sphere of influence where they can extract an- 
gular momentum from the binary [6lL |62(. A number 
of other mechanisms have been identified that can accel- 
erate the evolution of binary SBHs, even in collisionlcss 
nuclei [H3, Constraints on the binarity of the MW 
SBH are fairly tight |60| , and only one clear detection of 
a binary SBH has so far been made [Hll , suggesting that 
they may be rare. 

C. Dark matter 

Detectability of DM crests depends critically on the 
DM density at r^, roughly the outer boundary of the 
region in which the mass distribution is modified by the 
(single or binary) SBH. Traditionally there have been two 
approaches to estimating p x at the centers of galaxies; 
unfortunately they lead to rather different conclusions 
about p x (rh)- 

N-body simulations of gravitational clustering follow 
the growth of DM halos as they evolve via mergers in 
an expanding, cold-dark-matter (ACDM) universe. Halo 
density profiles in these simulations are well determined 
on scales 10~ 2 < r/r V i r < 10°, where the virial radius 
r V i r is of order 10 2 kpc for a galaxy like the Milky Way; 
hence inferences about p x on scales of w 10~ 3 kpc re- 
quire a radical extrapolation from the TV-body results. A 
standard parametrization of p x in these simulated halos 
is 

ft t (r) = Por 1 (l + 0" a (11) 

|66j . the "NFW profile," where £ = r/r s and r s is a scale 
length of order r v i r . In the Milky Way, r V i r ^> Rq (the 
radius of the Solar circle) hence Eq. l|llf) is essentially a 
power law at r < R@ and the implied DM density at 
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p x {r h ) « 3OM pc- 



10- 2 M o pc 



-3 



R, 



kpc 



3 pc 
(12) 

where p Q = p x (Rq) and p ~ 8 x 10 3 M pc 3 (from 
the Galactic rotation curve). Moore et al. argue 



for a steeper inner slope, p x 



-1.5 



implying p x {r h ) 



10 3 M Q pc~ 3 . Still higher DM densities could exist if the 
baryons (stars, gas) lose energy radiatively and contract, 
deepening the potential well and pulling in the DM (e.g. 
[E3)- Based on halo scaling relations (§3, DM densities 
at r/j would be higher in spheroids fainter than that of 
the Milky Way. 

Rotation-curve studies of dark-matter-dominated 
galaxies are generally inter pre ted as implying much 
lower, central DM densities [IS El El f 



While 

there are caveats to this interpretation - systematic bi- 
ases in long-slit observations |75| . non-circular motions 
[76| . gas pressure etc. - these effects do not seem 
capable of fully explaining the discrepancies between ro- 
tation curve data and expressions like l|ll(l |73l l78| . A 
model for p x (r) that is often fit to rotation curve data is 



p^^pcii+ny'ii+ey 



(13) 



the "Burkert profile,", where £ = r/r c and r c is the core 
radius. Inferred core radii are ~ 10 3 pc and inferred 
central densities typically lie in the range p c (1 — 5) x 
10- 2 M Q pc- 3 . 

Several resolutions have been suggested for this appar- 
ent conflict between theory and observation [j^. Since 
the TV-body halos are not resolved on the scales (~ 10 2 
pc) where rotation curves are typically measured, the 
mismatch may simply be due to a poor choice of empir- 
ical model used to parametrize p x (r). F° r instance, the 
alternative parametrization 



Px( r ) = P £ p exp -6£ 



(14) 



the "Prugniel-Simien" law 80], is both a better fit to 
the iV-body data than Eq. I|ll|) and is also in reasonable 
accord with rotation curve data [6j| 0, |s3i Ell • Here 
£ = r/R e and R e is the radius containing 1/2 of the 
projected halo mass (the relation between R e and the 
virial radius is discussed in HH). The profile shape in 
Eq. (|14fl is determined by the curvature parameter n; for 
the TV-body halos, n is found to be a weak function of 
halo mass, but exhibits a substantial scatter at all halo 



D. Summary of Galaxy Properties Relevant to the 
Formation of Crests 



Relaxation times are short enough for the formation 
of DM crests in stellar spheroids with L ra 1O 95 L or 
fainter. However the nuclear structure of faint galaxies 
is typically unresolved and relaxation times in spheroids 
with L < 1O 9 L are uncertain. SBHs appear to be ubiq- 
uitous in stellar spheroids brighter than ~ 1O 9 L , with 
masses that are well predicted by the properties of the 
stellar spheroid via empirical scaling relations. A handful 
of massive (~ 1O 5 M ) BHs have been detected in fainter 
galaxies via their non-thermal spectral features, but it is 
not clear what fraction of low-luminosity spheroids con- 
tain BHs. Dark matter densities at r ra r/j in spheroids 
with L ra lO 9 I/ may be as high as 10 2 — lO 3 Af pc 2 , 
if standard parametrizations of ACDM halo models are 
correct; or as low as ~ 10~ 2 — lO _1 M pc -3 if rotation 
curve studies are to be believed. 



V. OBSERVABILITY OF THE CRESTS 

In this section, we consider the implications of 
collisionally-generated DM crests for the rate of particle 
self-annihilations and for the detectability of the resul- 
tant gamma rays. We consider separately the case of the 
Milky Way and external galaxies. Since the p x ~ r -3 / 2 
crests considered here are relatively weak (e.g. compared 
with the steeper spikes that form via adiabatic growth of 
a SBH, p x ~ r~ 7 , 2 < 7 < 3 12]), we focus on the ques- 
tion of whether the DM distribution inferred above at 
r < implies a substantial increase in the predicted an- 
nihilation signal compared with the signal from a galaxy 
that lacks such a crest. 

The annihilation rate from neutralinos near the center 
of a spherically-symmetric DM halo is proportional to 
/ p x (r)r 2 dr. In the presence of a p x ~ r~ 3 / 2 crest, the 
integral diverges as log(r ( ^ 1 ) where ro is the inner radius 
of the crest. Roughly, ro is the maximum of (rs,r a ) 
where r$ = 2GM,/c 2 is the Schwarzschild radius of the 
SBH and r a is the radius where the self-annihilation time 
equals ~ T r . For all reasonable values of m x and av, 
the annihilation cross section, r a <C r$ hence we assume 
ro ~ rs in what follows. 

Evaluating this integral in the case of the TV-body den- 
sity profile of Fig. 3 at t = T r (rh), we find 



p 2 x (r)r 2 dr = p 2 h r\ 



f 1 In ( T ± ) C, 



(16) 



2 < n < 6 



(15) 



(e.g. |84j. Fig. la). The constants b and p in Eq. ljTI|) are 
determined uniquely by n [8j|; for the range in n cited 
above, 0.7 < p < 0.9, i.e. the density increases more 
slowly toward the center than the NFW [6(j (oc r -1 ) 
profile. 



with Ci w 4.22 x 10~ 2 , C 2 ~ 2.94 and p h = p x (r h ). 
(We assumed p x (r) = pi(r/ri)~ 3 / 2 at r < r\ = O.Olr/j 
where the A^-body density is poorly defined, as justified 
above.) Roughly 90% of the integral comes from matter 
at r <rh- 

While the DM density is affected by the SBH at all r < 
rn, the p x ~ r -3 / 2 crests only extend out to r < O.lr^. 
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We find that ~ 20% of the integral i|16|) comes from the 
crest proper, i.e. from r < O.lr^. 

Setting ln(rh/r ) w \n(c 2 / v 2 ms ) « 15, the integral in 
Eq. (|16fl is ~ 3.6p^r^. For comparison, the value corre- 
sponding to a constant density within is ~ 0.33p|r^. 

Henceforth we assume that Fig. 3 correctly represents 
the DM distribution at r < rh in nuclei that are ~ T r {rh) 
old, i.e. that 



p 2 x (r)r 2 dr^3p 2 h rl 



(17) 



The observable annihilation signal will typically include 
a contribution from r < , as well as a contribution from 
DM beyond that lies within the detector's window. We 
will consider several possible forms for the large-radius 
dependence of p x on r. 

At times later than ~ T r (r^), the amplitude of the 
crest drops due to continued heating from the stars (e.g. 
Fig. [21 • We discuss the consequences of this effect for the 
observability of crests in more detail below. 



A. Milky Way 

The photon flux from neutralino annihilations in the 
Galactic halo, observed by a detector with angular ac- 
ceptance Afl, is 



/(Am = — [ dn' [ P 2 Ji)di 
4?r j An j 



(18a) 



(18b) 



Here o~v is the annihilation cross section times relative ve- 
locity (in the nonrelativistic limit), dN/dE is the gamma 
ray spectrum per annihilation, and I is the line-of-sight 
distance. In the case of telescope centered on the Milky 
Way SBH, the integral can be written 



/(*) 



1 




R% 


/ Pl( r ) r 

- J r 






p x {r)r 2 dr 






fRo / 




p x [r -r 





* 2 i?l 



(19) 



where ^ 2 = AH/tt; the inequality reflects the omission 
of the contribution from DM outside the Solar circle, and 
we have also implicitly assumed ^i?© > r/,, valid for all 
current and planned detectors. (Atmospheric Cerenkov 
telescopes like HESS, and the proposed satellite observa- 
tory GLAST, have A!1k5x 10~ 5 sr.) We equate the 
first term on the RHS of Eq. JT5J} with Eq. JT7J. To eval- 
uate the additional terms we require p x {r) at r > r%. A 
simple model is a power-law, 




|0 910 r h/ R scal 



le 



FIG. 5: Contribution to the self-annihilation signal from DM 
at r < r^, i.e. within the SBH's influence radius. The DM 
distribution at r < th is assumed to be the same as in Fig. |3] 
at t — T r (rh)- Three different models for the DM density at 
r > Th have been evaluated: Lower solid line: Burkert profile; 
Upper solid line: NFW profile; Dotted lines: Prugniel-Simien 
profile, with n = (2,3,4,5,6), increasing upwards. R sca i e is 
the appropriate scale length for each model, i.e. r c (Burkert), 
r 3 (NFW), R e (Prugniel-Simien). 



For 7 = 1 this is a reasonable approximation to both 
an NFW profile, Eq. I|ll(l . and a Prugniel-Simien profile, 
Eq. ifHjl. at r < R®. Setting 7 = 1 gives 



/(*) 



P 2 x (r h )rlRo 2 



3 + - 



7T *i? Q 



(21) 



Px( r ) = Px( r h) (r/r h y 



r > r h . 



(20) 



where the first term is the contribution from r < 
(Fig. El. For Afl = 5 x 10~ 5 (* w 4 x 10~ 3 sr), roughly 
20% of I is due to matter at r < rj, and roughly 4% 
from the crest proper. This result suggests that the de- 
tectability of DM self-annihilations at the Galactic center 
is not likely to be significantly affected by the presence 
of a crest. We note that in Ref. ; 25j, the contribution of 
a DM crest at the Galactic center to the self-annihilation 
signal was found to be much greater. The discrepancy 
arises from the fact that these authors ignored the third 
(dominant) term on the right hand side of Eq. (|19|) . and 
they assumed that the crest extended as p x oc r -3 / 2 all 
the way to rh, rather than to only a fraction of rh- 

The relative contribution from the crest could be in- 
creased by assuming a steeper DM density falloff, e.g. in 
a halo where the DM had been pulled in by contraction 
of the baryons, or a smaller angular acceptance Af2 of 
the detector. As discussed above, observations favor a 
lower DM density implying even less of a contribution 
from the crest. 
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B. External Galaxies 

In the case of a detector centered on the SBH of an 
external galaxy, the geometrical term in the expression 
for the flux becomes 



1 

£2 



pi{r)r 2 dr 



(22) 



where D is the distance to the galaxy and the integral 
includes as much of the galaxy as is imaged by the tele- 
scope; thus 



1(9) 



1 



p 2 (r)r 2 dr + 



pt(r)r 2 dr 



(23) 



We again equate the first term in brackets with Eq. (|17|l . 
For the low-luminosity galaxies that are likely to harbor 
a DM crest (Sect. IVA), the gamma-ray telescope would 
image essentially the entire halo and so we set the upper 
limit of the second integral to infinity in what follows. 

Following the discussion in Sect. IIIC, we considered 
three possible forms for the DM density beyond r = r^: 
the NFW profile fEa.lTTj). the Burkert profile fEa. lT3|) 
and the Prugniel-Simicn profile fEci. I14|) . 

Fig. El shows the relative contribution of the DM at 
r < Th to the annihilation signal for the three assumed 
DM profiles. In the case of the Prugniel-Simien profile, 
the curvature parameter n has been varied over the range 
2 < n < 6 that approximately characterizes the ./V-body 
haloes ■ The contribution of the crest is insignificant 
in the case of the Burkert profile and only significant in 
the case of the NFW profile if r/j/r s is unphysically large. 
For the Prugniel-Simien profile however, the crest adds 
significantly to the total signal for large n and rh/R e > 
1CT 4 . 

Whether such large values of rh/R e are physically rea- 
sonable depends on the poorly-understood relations be- 
tween SBH mass, galaxy mass and DM halo properties at 
the low-mass end of these distributions. The scale radii 
i? e of A-body DM halos obey 



Re 



100 kpc 



1.1 



M, 



DM 



10 12 M G 



1/2 



(24) 



|69J. Thus 
rh_ ^ GNU 

Re ~ <T 2 R e 



4xl0- 6 ^ M ' 



M. 



DM 



V 10 6 M Q J V 100 km (TV V 1O 12 M 
Adopting the empirical relation between M, and cr* |To| 



M. 



5.7 



CT- k 



100 km s" 



a « 4.86 



(26) 



(which has only been established for er* > 100 km s 1 ) 
this becomes 



^«8x 10 
Re 



-«(. 



M. 



V f o 6 a/ q 



0.59 



M 



DM 



10 12 M a 



-1/2 



(27) 



If M, scales linearly with halo mass at the low-mass end 
of the distribution, this relation implies that rh/R e is 
essentially independent of M. and Mdm, hence roughly 
equal to its value in the Milky Way, r h /R e « I x 10 -5 . 

On the other hand, the empirical relation between ro- 
tation curve peak velocity and stellar velocity dispersion 
suggests a nonlinear relation between M, and Mdm, of 
the approximate form |85| 



M. 



10 6 M(t 



K 



M 



DM 



10 12 M ( 



with 1 .5 < (5 < 2 and 0.1 < K < I. Thus 



r h , ,0. 59-1/2/3 ,,0.59/3-1/2 
— (X M, (X M r 

Re 



DM 



(28) 



(29) 



and rh/Re would scale approximately as M, - 3 oc M^ M . 

These scaling relations should be considered highly un- 
certain for the reasons discussed in Sect. IV. A third pos- 
sibility, not inconsistent with the limited observational 
constraints on SBH and halo masses, is that SBH mass 
is essentially independent of halo mass at the low-mass 
end of the distributions. This assumption would imply 
increasing values of r^/Re in dwarf galaxies, hence larger 
contributions from the crests to the annihilation signal 
(Fig.©. 



C. Dependence on initial conditions and galaxy age 

The estimates made above of DM crest observability 
were based on the TV-body density profile, Fig. |3 at 
t = T r (rh). While this density profile is "universal" at 
r < O.Ir/j in the sense that p x ~ r~ 3 / 2 is a steady-state 
solution to the Fokker-Planck equation in a point-mass 
potential, the amplitude of the crest decays after its for- 
mation due to continued heating by the stars (Fig. 
[T7L |3). Furthermore the distribution of DM within 
the SBH's influence radius could be different if the initial 
conditions were very different from those assumed here. 

Initial conditions: Until now we have made the con- 
servative assumption that the DM density profile was 
initially very flat near the center. These are appropriate 
initial conditions if the galaxy hosting the DM experi- 
enced a merger following formation of SBHs, since a bi- 
^nary SBH is efficient at ejecting matter from the center 
of a galaxy and creating a low-density core 16]. However 
the mean time between mergers is a strong function of 
galaxy luminosity/mass, and an isolated, low-mass stellar 
spheroid might not have experienced a significant merger 
in the last I0 10 yr. (We note that this statement does 
not apply to the bulges of massive spiral galaxies, like 
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that of the Milky Way, since the merger probability is 
determined by the overall mass/radius of the galaxy, and 
these are much larger than the mass/radius of the bulge 
in the case of a spiral galaxy.) For instance, intermediate- 
mass black holes (IMBHs) might have formed via direct 
collapse of primoridial gas in low-mass halos |2£| ■ In this 
scenario, the distribution of DM around the IMBH could 
have the steep dependence with radius predicted by so- 
called "adiabatic growth" models, p x ~ r~ 7 ,2 < 7 < 3 
[87l HU . Such steep spikes would be "softened" by self- 
annihilations at small radii and by heating from the stars, 
but Fokker-Planck integrations show that the density 
profile after one relaxation time can remain considerably 
steeper than that of the collisional, p x ~ r~ 3 / 2 crests 
discussed here 0, El • The same considerations would 
apply to a spike in a more massive halo that did not hap- 
pen to experience a major merger since the epoch of SBH 
formation. 

Age: On time scales long compared with T r (r^), con- 
tinued collisional evolution of a nucleus can result in 
different stellar and/or DM density profiles at r < r/,. 
Heating of DM by stars causes the normalization of 
the DM density to drop, while roughly maintaining the 
p x ~ r -3 / 2 dependence at r <C r/, (e.g. Fig. 1, 2). Cap- 
ture or tidal disruption of stars by the SBH is effectively 
a heat source, causing a nucleus to expand |89j |: the time 
scale is of order ~ T r (r^) and depends on the sizes and 
masses of stars. If there is a range of stellar masses, the 
more massive stars will accumulate near the center, caus- 
ing the stellar mass density profile to become more cen- 
trally peaked [ll|. While the "steady-state" functional 
form for p x (r), Eq. (JJJ, is formally independent of p*(r), 
the rate at which the stars transfer energy to the DM 
depends on the stellar density. As discussed above (Sect. 
IV A), it is not clear that any nucleus harboring a SBH 
is much older than one relaxation time. 



VI. SUMMARY 

By considering the joint evolution of the stellar and 
dark-matter (DM) densities at the center of a galaxy con- 



taining a supermassive black hole (SBH), we have shown 
that the existence of a "crest" (collisionally regenerated 
dark-matter structure) can be robustly predicted in any 
nucleus old enough to have generated a Bahcall-Wolf cusp 
in the stars. This time scale is roughly 10 Gyr in the case 
of the Milky Way, and probably shorter for fainter galax- 
ies that contain massive BHs. Crest generation occurs 
even in the (probably generic) case of a nucleus that pre- 
viously experienced the scouring effect of a binary SBH. 
Standard galaxy scaling relations suggest that crests do 
not dramatically change the prospects for indirect detec- 
tion of DM. However this conclusion could be modified 
if the DM density falls steeply beyond the SBH's gravi- 
tational influence radius, or if SBH masses scale weakly 
with DM halo masses. In any case, crest formation im- 
plies a significant increase in the density of DM around 
the central SBH, and the presence of crests should be 
taken into account when studying the constraints on the 
annihilation rate from stellar orbits 90], or the evolu- 
tion of stars in the innermost regions of galactic nuclei 

MM 
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